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Abstract 

This paper summarizes some work I've been doing on eigenvalue 
correlators of Random Matrix Models which show some interesting 
behaviour. First we consider matrix models with gaps in there spec- 
trum or density of eigenvalues. The density-density correlators of these 
models depend on whether N, where N is the size of the matrix, takes 
even or odd values. The fact that this dependence persists in the 
large N thermodynamic limit is an unusual property and may have 
consequences in the study of one electron effects in mesoscopic sys- 
tems. Secondly, we study the parametric and cross correlators of the 
Harish Chandra-Itzykson-Zuber matrix model. The analytic expres- 
sions determine how the correlators change as a parameter (e.g. the 
strength of a perturbation in the hamiltonian of the chaotic system or 
external magnetic field on a sample of material) is varied. The results 
are relevant for the conductance fluctuations in disordered mesoscopic 
systems. 
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1 Introduction 



Recently there is a lot of activity in the field of quantum chaos and mesoscopic 
systems. This has resulted in the study of eigenvalue correlators in large ran- 



dom matrix models ref. p|-flO|. This paper summarizes some work I've been 



doing on calculating eigenvalue correlators in random matrix models which 
show novel properties which maybe useful in the study of transport proper- 
ties in mesoscopic conductors. The first models are random matrix models 
with gaps in the eigenvalue spectrum and the second are two-matrix mod- 



els first considered by Harish Chandra-Itzykson-Zuber ref. |TT]. In the first 
model the 'fine grained' correlators are found using the method of orthogo- 
nal polynomials ref.[|13|]. These correlators are unusual in that in the large N 
thermodynamic limit they tend to different limits depending on whether N 



goes to infinity through even or odd ref. [Q. This property maybe found in 
mesoscopic systems which are sensitive to single electron effects. The second 
models are the two-matrix models of Harish Chandra- Itzykson-Zuber ref. 
||, |J. Here correlators are calculated using the Dyson- Schwinger method 
and give the long ranged parametric and cross-correlators. Transport exper- 
iments involving changes in magnetic fields in which long range eigenvalue 
statistics are effected will be a testing ground for these results. 

I shall first discuss the double well matrix model with potential V(M) = 
— ^M 2 + ^M 4 where M is a random NxN matrix. I will present some results 
for the two-point density- density correlation functions which show interest- 
ing parity effects and further characterizes these models in a new universality 
class. The original model is Z2 symmetric while there is Z2 symmetry break- 
ing in the correlation functions. Second, the Chandra-Itzykson-Zuber matrix 
model is discussed and the results for the smoothed long range parametric 
and crossover eigenvalue correlators, found using the Dyson-Schwinger equa- 
tions, are given. 



2 Notations and Conventions for the Double- 
Well Matrix Model 

We start by establishing the notations and conventions. Let M be a hermitian 
matrix. The partition function to be considered is Z = / dMe~ NtrV ^ M ^ where 
M = NxN hermitian matrix. The Haar measure dM = rjJLi dM u Ui<j dM^fdM^ 
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with My = M\j + iM-j and A^ 2 independent variables. V(M) is a poly- 
nomial in M: V(M) = 9l M + (g 2 /2)M 2 + (g 3 /3)M 3 + (^ 4 /4)M 4 + The 

partition function is invariant under the change of variable M' = UMU^ 
where U is a unitary matrix. We can use this invariance and go to the 
diagonal basis ie D' = UMU^ such that D' is the matrix diagonal to M 

with eigenvalues Ai,A2, Xn- Then the partition function becomes Z = 

C J" oo n^irfA J A(A) 2 e- 7V ^-i v ' (A!) where A(A) = Ui<j W ~ A; I is the Van- 
dermonde determinant. The integration over the group U with the appro- 
priate measure is trivial and is just the constant C. By exponentiating the 
determinant as a 'trace log' we arrive at the Dyson Gas or Coulomb Gas 
picture. The partition function is simply Z = C /f^ n£Li d\e~ s ^ with 
S(X) = jV£^ V(Xi) - 2£ M .^HA l - A,|. 

This is just a system of N particles with coordinates Aj on the real line, 
confined by a potential and repelling each other with a logarithmic repulsion. 
The spectrum or the density of eigenvalues p(x) = j? J2iLi ${ x ~ Aj) is in the 
large N limit or doing the saddle point analysis just the Wigner semi-circle for 
a (Gaussian probability distribution for the eigenvalues) quadratic potential. 
The physical picture is that the eigenvalues try to be at the bottom of the 
well. But it costs energy to sit on top of each other because of logarithmic 
repulsion, so they spread, p has support on a finite line segment. This 
continues to be true whether the potential is quadratic or a more general 
polynomial and only depends on there being a single well though the shape 
of the Wigner semi-circle is correspondingly modified. For the quadratic 
potential the density is p(x) = — a){b — x) where [a, b] are the end of 

the cuts. See Fig. [I]. 

On changing the potential more drastically by having two humps or wells, 
the simplest example being a potential V(M) = — ^M 2 + |M 4 , the density 
can get disconnected support. The precise expressions for the density of 
eigenvalues are as follows: 



p{x) = —xJ(x 2 — a?){b 2 — x 2 ) a<x<b 

= -b<x<-a (2.1) 

where a 2 = -[\fJ>\ — and b 2 = -[|//| + 2 v /g] and \p\ > 2 v /g, which 
is the condition that the wells are sufficiently deep. The eigenvalues sit in 
symmetric bands centered around each well. Thus p has support on two line 
segments. As \p\ approaches 2y/g, a — > and the two bands merge at the 
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Figure 1: (a). The confining potential (b). The density of eigenvalues 



origin. The density is then 



gx A 2fi 
p(x) = \ x z 

7T V 9 



I2\n\ 



< x < 



I2\n\ 







otherwise. 



(2.2) 



The phase diagram and density of eigenvalues for the M 4 potential is 
shown in Figs. ^| 

The simplest way to determine p(z) explicitly is to use the generating 
function F(z) =< jfTr-^jj > and its saddle point or Schwinger-Dyson 
equation also known in the mathematics literature as the Riemann-Hilbert 
problem F(z) = \[V\z) + VA] with A(z) = V'(z) 2 - Ab(z) and b(z) = 
gz 2 + fj, + g < -^TrM 2 > (see ref. [0]). The density p{x) is then determined 



by the formula p(z) = —^lm^jA(z). 

In what follows I will outline the recurrence coefficient method of the 
orthogonal polynomials for the two-cut matrix model see ref. [12| for more 
details. I give the results for the two-point correlators (also known as the 
'smoothed' or 'long range' correlators) for all three ensembles. Then the two- 
matrix model of Chandra-Itzykson-Zuber is defined and the expressions for 
the parametric and crossover correlators are given. 
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Figure 2: (a). The double- well potential (b). Density of eigenvalues (c). 
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3 Orthogonal Polynomial Approach 



The partition function Z, can be rewritten in terms of the orthogonal polyno- 
mials P n where the polynomials are defined as /f^ d\e~ NV ^ x > P n (X)P m (X) = 
hj nm then P„ = A n + l.o. and P (X) = 1, Pi (A) = A + c, P 2 (A) = A 2 + aA + 6 

Then Z the partition function in terms of the orthogonal polynomials 

is 

P (Ai) Po(A^v) 
Pi(Ai) .... Pi(Aiv) 



N 



-NJ2V(X) 



1=1 



(3-3) 



Pjv-i(Ai) Pzv-i(Ajv) 

The partition function Z is also known if we know the tin's as the partition 
function can be expressed in terms of the h n 's Z = N\h hih 2 h^-i- For ex- 
ample: N=2 case / dX 1 dX 2 e- NV ^- NV ^(P (X 1 )P 1 (X 2 ) - P (A 2 )P 1 (A 1 )) 2 = 
h hi + h hi = 2!/i /ii- So the question is how does one find the /i's? 
The P n satisfy recurrence relations 

(3.4) 

Note that J xP n P n ^ 2 e~ NVix) dx = as xP n _ 2 = P n _i + l.o.. Thus P n _ 2 & 
l.o. terms do not appear on the right hand side of the recurrence relation eq. 
(|3.4|). Then as h n = h n -iR n the product h hi fo/v-i — h (h R 1 )(h RiR2) 



.(h Ri...R N ) = h^R^- l R^-' A Ptv-i- The free energy T = InZ 



IniV! + iVln/io + J2n=i(N ~ n ) ^-Rn hence we need recurrence coefficients 
P n 's to get the free energy. 



Asymptotic Ansatz for the Orthogonal Poly- 
nomials of the Symmetric Double- Well Ma- 
trix Model 



We have been able to derive in ref. |13| for the symmetric double-well ma- 
trix model the asymptotic ansatz for the orthogonal polynomials ip n {X) = 



^-exp(-fV(X)) which is 



^n(A) 



/(A) 



[cos(iVC - {N - n)<j> + x + (- 



■l)^)(A)+0(-)] (4.5) 
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where /, (, 0, x an d T] are functions of A, and ip n is damped outside of the cuts. 
We show that /(A) = ^ ^ ~ a sin 20(A) from the orthogonality condition 
satisfied by the orthogonal polynomials. C'{X) = —7r P(^) from the relation 
K(\, A) = p(A) where K(/j,u) is the kernel defined by ^ X^I 1 '^(/-O'fM^) 
and determines all eigenvalue correlators. While 0(A) and 77(A) are deter- 
mined from the recurrence relations satisfied by the orthogonal polynomials 

cos20(A) = X2 pCjf\ cos2r ? (A) = and sin2r ? (A) = 

/, 0, 77 and x are universal functions independent of the potential, the 
only dependence on V enters through the endpoints of the cuts a and b. ( is 
non-universal since the eigenvalue density depends in general on the detailed 
form of V. More details are to be found in the ref. fl3], [I4l| . From this one 
can establish that the gapped matrix model is in a new universality class. 

Using the asymptotic ansatz the full density-density correlation function 
maybe found and is given in ref. [TB[]. A simpler result are the 'smoothed' or 



'long range' two-point correlation functions found in ref. |14| using a method 
of steepest decent which has a constant C which is undetermined unless one 
is explicitly in the symmetric case where C = (— 1)^ ref. (ref. |L6| obtain 
this result using a method due to Shohat), other values for C have been found 
earlier using the loop equation method ref. ||15|| . Recently in ref. |T7| it was 
shown that for double wells with equal depths but unequal widths, in the 
limit of the symmetric double wells, give the value C = (— 1)^ confirming 
the results of ref. 



ln 2 A 2 //4A./0 - 



{ (A-^ + (A-„) +A +/i -2 {X + fi) + 2C )- (46) 

Here a(z) = (z 2 —a 2 )(z 2 — 6 2 ), s = ai+a 2 +a 3 +a^, e\ = +1 for a 3 < A < 
^4, £a = —1 for ai < A < a 2 and j3 = 1,2,4 depending on whether M the 
matrix is real orthogonal, hermitian or self-dual quartonian. This result is 
different for even and odd N and hence has broken the Z 2 symmetry. It would 
be very interesting to be able to see this effect in experiments of mesoscopic 
systems which are sensitive to single electron effects. 
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5 The Harish Chandra-Itzykson-Zuber ma- 
trix model and density-density correlators 

The two-matrix model of Harish Chandra-Itzykson-Zuber is defined by the 
partition function 

Z = J dAdBe- s (5.7) 

where S = NSp[V(A) + V(B) - cAB] and V{A) = \pA 2 . For A = ±, 
A and B are N x N real symmetric matrices (orthogonal ensemble), for 
A = 1, A and B are N x N Hermitian matrices (unitary ensemble), and 
for A = 2 N x N real self-dual quaternions (symplectic ensemble). Sp(A) 
stands for TrA for A = ~, 1 and for \TrA for A = 2. In this model the 
connected density-density correlator is pab(x,u) =< Pa(x)pb(u) >c, where 
the density is defined as Pa( x ) = jjTr8(x — A), < X >= A / dAdBe~ s X , 
and the subscript c implies connected part. In ref. [§] the Dyson- Schwinger 
equations are used to derive these eigenvalue correlators. Here only the 
results and physical interpretations are elaborated on. 

In the large-N limit, the expectation value for the density is given by the 
well known Wigner semicircle law 



< Pa(x) >=< Pb{%) >= — o v/fl2 x2 ' l x l — a ' (^-^) 



ira 2 



and < p >= for \x\ > a, where a, the "end point of the cut" is given 
by a = ( The result for the connected density-density correlator to 

leading order in valid over the entire cut, is 



. Ill 

PAB{x,y) = — x 

4tt z N z Xa z cos cos a 

1 + cosh u cos(9 + a) 1 — cosh u cos(9 — a) 



[cosh u + cos(^ + a)] 2 [cosh m — cos(6 ) — a)] 2 ' 



(5.9) 



where u = ln(— ), |y| < a, and we have defined sin 6 = - and sin a = K . 
For A = 1 the above result was derived in ref. using different methods 
from the Dyson-Schwinger method. The Dyson- Schwinger method is capable 
of generalization to A = ~,2 with the above result. Then these results are 
relevant for the calculation of conductance fluctuation of mesoscopic systems 
in which the magnetic field is changing. 
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When one is interested in transitions from one symmetry class to another, 
a hamiltonian A is considered consisting of two parts B and V each drawn 
from different ensembles A — B + V. The partition function Z and action S 
are modified to have 

V(A) = ^A 2 , V(B) = ^ 2 B 2 (5.10) 

and 

fJ.ijl-fJ.iN) 

c = Mi, Z*2 = (2 _ ^ ■ (5-11) 

The constant /ii measures the strength of the perturbation. At \i\ = oo 
we get GOE, while A*i = ^ gives GUE. However, we will work in the more 
general case where fix, \±i and c are independent 0(1) parameters. This is 
just the standard two-matrix model except that A and B are drawn from 
different ensembles. The choice of parameters eq. ( 5.11|) , mentioned above 



corresponds to the crossover from GOE to GUE is then a special case of the 
formula derived. 

We are interested in calculating the connected density-density correlator 

p c aa( x , y) =< pA(x)p A (y) > c (5.12) 

The full smoothed global result for the connected eigenvalue correlator is (see 
ref. i) 

p c aa( x i y) = paa( x > y) + paa*( x i y) ( 5 -i3) 

where we find 

I ( \ 1 1 (a 2 — xy) , . 



and 



ii , v c 2 /i 2 r cos(0 - 9) - 1 



N 2 (2hi/j, 2 — c 2 ) 2 cos 4> cos 9 

2/i 2 (x - y) 2 - Jr (cos 20 + cos 2# - 2 cos(0 + 0)) - ^(x - y)(sin0 - sin0) 

(/* - a yf + J(l - cos(£ - 0)) - ^(x - y)(sm6 - sin0) 2 ) 

cos(0 + 0) + l 
H x 

2 cos (f) cos 

2/i 2 (x-y) 2 - ^(cos20 + cos2g + 2cos(g + 0)) - ^(x - y) (sing - sin 

(fJ 2 i(x -yf + J(l + cos(# + 0)) - if (x - yKsine - sin0) 2 ) ^' 15) 



where we have used sin 9 = ~, sin <b 
notice that for \i\ = oo, 



y - and a 2 = ±. 

a a 



After some algebra we 



Paa( x ,v) 



-1 



(5.16) 



7r 2 N 2 (x — y) 2 



which is the GOE result and for fii 



2_ 
N ■ 



p c aa(z, y) 



-l 



(5.17) 



2n 2 N 2 (x -y) 2 



the GUE result. The expression eq. ( |5.15|) is relevant for crossover from 
the GOE to GUE ensemble. One application of these correlation functions 
is to disordered mesoscopic systems using the transmission matrix formalism 
and the other is in the study of unoriented random surfaces. 

6 Conclusions 

In conclusion we have presented two classes of random matrix models. One 
in which there are gaps in the eigenvalue distribution and the other in which 
there are two coupled matrices drawn from the three ensembles (matrices 
taken from the same and different ensemble have been considered). In each 
of the models we have derived eigenvalue correlators particularly density- 
density correlators. In the first case of gapped matrix models we have eigen- 
value correlators which dependent on N the size of the matrix. This be- 
haviour persists in the large N thermodynamic limit and for the symmetric 
double- well matrix model parity effects are present. For the coupled matrix 
models long range smoothed correlators are found. These are the paramet- 
ric and crossover correlators which maybe found in mesoscopic experiments. 
Density-density correlators are applicable in calculations of conductance fluc- 
tuations of mesoscopic conductors. Our results in these models are valid for 
all eigenvalues near the center as well as the edge of the semi-circle. The be- 
havior near the edge of the cut is particularly relevant in studies of transport 
properties of mesoscopic conductors ref. p. Thus clever mesoscopic exper- 
iments should be devised which will show the effects found in both types of 
these matrix models. 
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